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Abstract 

We show that, in spite of a rather common opinion, quantum me- 
chanics can be represented as an approximation of classical statistical 
mechanics. The approximation under consideration is based on the 
ordinary Taylor expansion of physical variables. The quantum contri- 
bution is given by the term of the second order. To escape technical 
difficulties related to the infinite dimension of phase space for quan- 
tum mechanics, we start with a detailed presentation of our approach 
for the finite dimensional quantum mechanics. We also separate real 
and complex cases, because the reproduction of the complex structure 
of quantum mechanics is a special problem which is not related to ap- 
proximation of classical averages. In our approach quantum mechanics 
is an approximative theory. It predicts statistical averages only with 
some precision. In principle, there might be found deviations of av- 
erages calculated within the quantum formalism from experimental 
averages (which are supposed to be equal to classical averages given 
by our model). 

Keywords: quantum and classical averages, von Neumann trace for- 
mula, approximation, small parameter, Taylor expansion 
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1 Introduction 



The problem of coupling the quantum probabilistic model (which is based 
on the Hilbert space calculus) and the classical probabilistic model (which is 
based on the measure-theoretic calculus) was intensively discussed already by 
fathers of quantum mechanics, see, e.g., the correspondence between Einstein 
and Schrodinger jlj. This is the problem of huge complexity, see, e.g., [2] [H] 
for different viewpoints and debates. 

Now days there is a rather common opinion that the probabilistic struc- 
ture of quantum mechanics cannot be considered simply as a special math- 
ematical representation of classical (measure-theoretic) probability theory. 
Such an opinion is based merely on a number of no-go theorems, see appendix 
2 for details. On the other hand, there are known various prequantum models 
that reproduce (at least some) features of quantum mechanics: De Broglie's 
double solution theory, Bohmian mechanics, Nelson's stochastic mechanics, 
t' Hooft's deterministic discrete models, see jB] -JD]- All such models are 
either nonlocal (as Bohmian mechanics and Nelson's stochastic mechanics) 
or reproduce only some (not all) predictions of quantum mechanics. In any 
event they do not contradict no-go theorems. 

In author's papers [TTJ , [TJ] there was proposed a new prequantum model: 
Prequantum Classical Statistical Field Theory, PCSFT. This model is very 
natural, because this is nothing else than standard classical statistical me- 
chanics: a) state space is phase space; b) variables are functions on phase 
space; c) statistical states (describing ensembles of systems) are probability 
measures. One important point is that phase space is infinite-dimensional. 
Prequantum states ("hidden variables") can be represented as vector fields, 
ijj(x) = (q(x),p(x)). However, this is merely a technical mathematical feature 
of PCSFT. The crucial point is a tricky way in which PCSFT is projected 
onto quantum mechanics. This projection is asymptotic. Quantum mechan- 
ics can be considered as an approximative theory: the quantum average given 
by the von Neumann trace formula appears as the first nontrivial contribu- 
tion into the classical average (given by the Lebesgue integral). 

To clarify the main distinguishing features of our theory, PCSFT, we shall 
divide its presentation into a few steps. First we consider the finite dimen- 
sional case. Here all computations are reduced to asymptotic expansions of 
simple Gaussian integrals. Then we shall proceed to the infinite-dimensional 
case where we should consider Gaussian integrals over functional spaces. We 
also start with consideration of a toy-model of quantum mechanics over reals, 
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see also [TTj, [T2]. This will simplify the understanding of PCSFT, because 
the prequantum phase-space model inducing the quantum model over com- 
plex numbers has a rather nontrivial geometric structure. And its presenta- 
tion should be separated from asymptotic expansion of Gaussian integrals. 

We remark that construction of a prequantum classical statistical model 
for finite-dimensional quantum mechanics is interesting not only from purely 
mathematical viewpoint. Since electron-spin and photon-polarization can be 
described by two dimensional complex spaces, in the finite-dimensional case 
our approach shows that that its is possible to construct a pure classical 
phase-space models for these quantum phenomena, see appendix 1. 

In our approach quantum mechanics is an approximative theory. It pre- 
dicts statistical averages only with some precision. In principle, there might 
be found deviations of averages calculated within the quantum formalism 
from experimental averages (which are supposed to be equal to classical av- 
erages given by our model). But at the moment our predictions is not of a 
high value for experimentalists, because PCSFT does not predict the mag- 
nitude of a small parameter a in the asymptotic representation of classical 
averages. In [TT], [T2] we speculated that the small parameter of PCSFT 
a (the dispersion of prequantum fluctuations) can be chosen equal to the 
Planck constant h. But that speculation was not justified. We notice that 
our asymptotic considerations are totally different from the standard con- 
siderations on the classical limit of quantum mechanics: obtaining classical 
phase space mechanics as the limit of quantum mechanics when the Planck 
constant h (which is considered as a small parameter) goes to zero. In our 
approach when the small parameter a goes to zero we obtain quantum theory 
as the limit case of classical (and not vice versa). In particular, neglecting 
by h induces neglecting by "quantum rotation", spin. But in PCSFT such 
degrees of freedom are not neglected, see appendix 1 for details. 

2 The Taylor approximation of averages for 
functions of random variables 

Here we follow chapter 11 of the book of Elena Ventzel. This book was 
written in the form of precise instructions which student should follow to 
solve a problem: 

"In practice we have very often situations in that, although investigated 
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function of random arguments is not strictly linear, but it differs practically 
so negligibly from a linear function that it can be approximately considered as 
linear. This is a consequence of the fact that in many problems fluctuations of 
random variables play the role of small deviations from the basic law. Since 
such deviations are relatively small, functions which are not linear in the 
whole range of variation of their arguments are almost linear in a restricted 
range of their random changes," ^3], P- 238. 

Let y = f(x). Here in general / is not linear, but it does not differ so 
much from linear on some interval [m x — 5,m x + S], where x = x(oj) is a 
random variable and 

m T = E x = I x(u) dPiuj 



is its average. Here 5 > is sufficiently small. Student of a military college 
should approximate / by using the first order Taylor expansion at the point 
m x : 

y(u) « /K) + f{m x )(x{u)) - m x ). (1) 
By taking the average of both sides he obtains: 

m y « f{m x ). (2) 

The crucial point is that the linear term f'(m x )(x(u) — m x ) does not give 
any contribution! Further Elena Ventzel pointed out [13 , p. 245: "For some 
problems the above linearization procedure may be unjustified, because the 
method of linearization may be not produce a sufficiently good approxima- 
tion. In such cases to test the applicability of the linearization method and 
to improve results there can be applied the method which is based on pre- 
serving not only the linear term in the expansion of function, but also some 
terms of higher orders." 

Let y = f(x). Student now should preserve the first three terms in the 
expansion of / into the Taylor series at the point m x : 

y(w) w f(m x ) + f'(m x )(x(uj) - m x ) + ~f"(m x )(x(uj) - m x f . (3) 



Hence 



m y sa f{m x ) + -£f"{m x ), (4) 
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where 



a 2 = E (x — m x ) 2 = J (x(uj) — m x ) 2 dP(u) 

is the dispersion of the random variable x. 

Let us now consider the special case of symmetric fluctuations: 



m x = 



and let us restrict considerations to functions / such that 



/(0) = 0. 



Then we obtain the following special form of (@J): 



m 



y 



f/"(0). 



(5) 



We emphasize again that the first derivative does not give any contribution 
into the average. 

Thus at the some level of approximation we can calculate averages not by 
using the Lebesgue integral (as we do in classical probability theory), but by 
finding the second derivative. Such a "calculus of probability" would match 
well with experiment. I hope that reader has already found analogy with the 
quantum calculus of probabilities. But for a better expression of this analogy 
we shall consider the multi-dimensional case. Let now x = (x\, ...,x n ), so we 
consider a system of n random variables. We consider the vector average: 
m x = (m Xl , ...,m Xn ) and the covariance matrix: 



By using the Taylor expansion we would like to obtain an algorithm for 
approximation of the average m y . We start directly from the second order 
Taylor expansion: 



b x = (irj). /;;:' /• (,-, 



m Xi ) (xj - m x .). 



We now consider the random variable 



x n {uj)). 



y(u) w /K, ...,m : 




m Xl ,...,m Xn )(xi(u) - m Xi ) 
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1 d 2 f 

+ 9 Yl F) F) i , m Xl ,...,m Xn )(x i (uj)-m Xi )(x j (uj)-m Xj ), (6) 



2 ^ dxidxj 

1,3=1 ■> 

and hence: 

1 n d 2 f 

m y « f(m xi ,...,m xi ) + ^Yl ft^T^ 1 ' ""' m *^ B x ' ( 7 ) 

i,j=l 1 3 

By using the vector notations we can rewrite the previous formulas as: 

y(uj) « /(m x ) + (/'(m x ),x(cj)-m x ) + i(/ // (m :E )(a;(cj)-m :E ),a;(cj)-m :E ). (8) 
and 

m» ~ /W + ^Tr B x f"(m x ). (9) 
Let us again consider the special case: m x = and /(0) = 0. We have: 

m y « ^Tr B x f"(0). (10) 

We now remark that the Hessian f"(0) is always a symmetric operator. Let 
us now represent / by its second derivative at zero: 

f^A = f"(0). 

Then we see that, at some level of approximation, instead of operation with 
Lebesgue integrals, one can use linear algebra: 

my « ^Tr B X A (11) 

We now proceed in mathematically rigorous way, namely, we shall es- 
timate the reminder which was neglected in the approximative formula for 
average. We also formalize correspondence between classical and quantum 
statistical models. 
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3 Classical and quantum statistical models 



3.1 Classical statistical model 

Classical statistical mechanics on phase space ^2?! — R n x R ri can be consid- 
ered as a special classical statistical model. In general a classical statistical 
model is defined in the following way: 

a) . Physical states uj are represented by points of some set Q (state space). 

b) . Physical variables are represented by functions / : Q — > R belonging 
to some functional space VlQ). 1 

c) . Statistical states are represented by probability measures on Q be- 
longing to some class S(Q). 

d) . The average of a physical variable (which is represented by a function 
/ G V(Q)) with respect to a statistical state (which is represented by a 
probability measure p G S(Q)) is given by 



In classical statistical mechanics Q = Vt 2n is phase space, V{Vt2 n ) = 
2n) is the space of all smooth functions on phase space, S(D,2 n ) is the 
space PM(f2 2 „) of all probability measures on phase space and the average 
is given by the Lebesgue integral on the a-algebra of Borel subsets of Q^n- 

Remark 3.1. We emphasize that the space of variables V(Q) need not 
coincide with the space of all random variables RV(Q) - measurable functions 
£ : f2 — > R. For example, if £1 is a differentiate manifold, it is natural to choose 
V(fi) consisting of smooth functions; if is an analytic manifold, it is natural to 
choose V(fi) consisting of analytic functions and so on. The space of statistical 
states S(Q) need not coincide with the space of all probability measures PM{Q). 
For example, for some statistical model S(fl) may consist of Gaussian measures. 

1 The choice of a concrete functional space V(Cl) depends on various physical and math- 
ematical factors. 




(12) 



A classical statistical model is a pair 



m = (s(n),v(n)). 
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3.2 Real finite-dimensional quantum mechanics 

We shall use a toy model of quantum mechanics which based on the real 
space. Statistical features of the correspondence between a prequantum clas- 
sical statistical model and quantum mechanics are more evident for this toy 
model. Denote the algebra of all (m x m) real matrices by the symbol 
M^ r \m). We denote by P"(m) the class of nonnegative symmetric trace- 
one matrices p G M^ r \m). We call them "density operators." We denote by 
L^\m) the class of all symmetric matrices. In the quantum model (for the 
m-dimensional real space) statistical states (describing ensembles of systems 
prepared for measurement) are represented by density matrices and quantum 
observables by matrices belonging Cs (m). The quantum average of an ob- 
servable A G Cs\m) with respect to a statistical state p G V^(m) is given 
by the von Neumann trace class formula |3]: 

< A > p = Tr pA. (13) 

In the operator representation observables and density matrices are corre- 
sponding classes of R-linear operators. Denote the quantum model by 

If m = 1, then quantum observables are given by real numbers (operators 
of multiplication by real numbers on the real line) and there is only one 
statistical state p = 1. Here < A > p = pA = A. 

3.3 Complex finite-dimensional quantum mechanics 

Denote the algebra of all (m x m) complex matrices by the symbol M^ c \m). 
We denote by T>^ c \m) the class of nonnegative symmetric trace-one matrices 
p G M^(m). We call them "density operators." We denote by d c \m) the 
class of all symmetric matrices. In the quantum model (for the m-dimensional 
complex space) statistical states (describing ensembles of systems prepared 
for measurement) are represented by density matrices and quantum observ- 
ables by matrices belonging ci c \m). The quantum average is given [3] by 
(|T3*|) . In the operator representation observables and density matrices are 
corresponding classes of C-linear operators. Denote the quantum model by 

< ) ant = (^ (C) M,4 C) M)- 
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If m = 1, then quantum observables are given by real numbers (operators of 
multiplication by real numbers on the complex plane) and and there is only 
one statistical state p — 1. Here < A > p = pA = A. 

3.4 Complex quantum mechanics 

Denote by H c a complex (separable and infinite-dimensional) Hilbert space. 
Denote the algebra of all bounded operators A : H c — > H c by the symbol 
C(H C ). The real linear subspace of C(H C ) consisting of self- adjoint operators 
is denoted by the symbol C S (H C ). We denote by T>(H C ) the class of non- 
negative trace-one operators p G C S (H C ). These are von Neumann density 
operators. In the quantum model statistical states (describing ensembles 2 
of systems prepared for measurement) are represented by density operators 
and quantum observables by operators from £. S (H C ). The quantum average 
is given jlj by (JT3J). Denote the quantum model by 



4 Taylor approximation of classical averages: 
one dimensional case 

States of systems are represented by real numbers, q G Q = R. Ensembles of 
such systems are described by probability measures on the real line, statistical 
states. We consider a special class of preparation procedures. They produce 
ensembles of systems described by Gaussian probability distributions on Q 
having the zero mean value and dispersion 



where as always |0(a 2 )| < Ca 2 for some constant C and a sufficiently small 
a. The crucial point is that a is a small parameter of our model. Denote this 
class of probability distributions by the symbol Sq(Q). 

2 We follow so called ensemble interpretation of quantum mechanics: Einstein, Mar- 
genau, Ballentine, Balian, Nieuwenhuizen and many others, see, e.g., |14| . By such an 
interpretation even a pure state (normalized vector of H c ) represents an ensemble. In 
the orthodox Copenhagen interpretation a pure state represents the state of an individual 
system, e.g., electron. 



iV quant (# c ) = (V(H C ),£ S (H C )). 




(14) 
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For a probability /x G Sq(Q), we have: 



, 2 



e 2(a + 0(c^)) ^qt 

d/i(g) = =. (15) 

v/27r(a + 0(a 2 )) 

We recall that, for a probability with the zero mean value, its dispersion is 
given by 

1 f°° -i 2 

= Aw r^7 2 V T / q 2 e^+°^dq. (16) 

As was already pointed out, we consider a as a small parameter. Therefore 
Gaussian probability distributions are very sharply concentrated around the 
point q = 0. By using the terminology of functional analysis we say that {/i = 
/i(a)} is a 5-family: lim a ^ M a ) = <5 in the sense of theory of distributions. 

In the approximation a = all systems are located at a single point, 
namely, go- However, a finer description (in that a can not be neglected) 
provides the picture of Gaussian bells concentrated nearby go- We remark 
that in average a system cannot go far away from go. By using the Chebyshov 
inequality one obtain for any C > : 

li{q : |g| >C}< a + ^ 2) - O, a - 0. (17) 

But the probabilistic inequality (|17|) does not exclude the possibility that 
some system could move far from go (of course, with a small probability). 

We also introduce a class of physical variables in the classical statistical 
model under consideration: 

a) / G C°°(R), a smooth function; 

b) /(0)=0; 

c) |/W(g)| < C/ e r /M c/,r/>0. 

Denote this functional space by the symbol V(Q), Q — R. 

The restriction to the growth of the fourth derivative will be used when 
we shall consider the Taylor expansion of / up two the fourth term. The ex- 
ponential growth implies integrability with respect to any Gaussian measure. 

Lemma 4.1. Let f G C n (so it is n times continuously differentiable) 
and let its nth derivative has the exponential growth. Then all derivatives of 
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orders n = 0, ...,n — 1, also have the exponential growth (in particular, f(q) 
grows exponentially). 

Proof. Under these conditions we can use the Taylor expansion with the 
integral remainder: 

/(?) = 7(0) + f(0)q + f -^- + + 9)^f^\9q)dB. 

(18) 

Since the growth of any polynomial can be compensated by decreasing of the 
e~ r M, by using the exponential estimate for the nth derivative we obtain: 

\f(q)\ = C ie rM + cS, [ (1 - er^e^de < Ce r| < (19) 
n ' Jo 

Here all constants depend on /. 

This simple exercise from the course of analysis will be useful in our 
further considerations. We defined the following classical statistical model 
on the real line: 

A) . States of systems are real numbers. 

B) . Statistical states (ensembles of systems) are represented by Gaussian 
probabilities having zero average and dispersion <r 2 (u) = a + 0(a 2 ), a — > 0. 

C) . Physical variables are smooth functions with exponentially growing 
fourth derivative which map zero into itself. 

We denote this model by A^ ass = (Sq(Q),V(Q)). As always in classical 
statistical physics, the average of a physical variable / e V(Q) with respect 
to an ensemble of systems which is described by a probability u e Sq(Q) is 
given by the integral: 

< f To ( , nf 2^ / fi^^^dq. (20) 

Since a is a parameter of the model, we can consider averages as functions 
of a :< / > M =< / > M (a). We are interested in the asymptotic expansion of 
averages when a — > 0. In particular, such an asymptotic expansion will give 
us the possibility to calculate averages approximately. 

Lemma 4.2. Let f e V(Q) and let u G 5g(Q). T/ien 

</>,(«) = f/"(0) + O(a 2 ). (21) 
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Proof. We start with the scaling of the state variable: 

q = a(n)x (22) 

We have: 

< / >, (a) = -= = / f(a(rix)e=t-dx. (23) 

V27T J-oo 

We now expand f(a(ii)x) by using the fourth order Taylor formula with the 
integral remainder, see Lemma 4.1: 

< / >„ («) = (0) (24) 



+|^| J^j\j\l-9f^\a^)x9)d9)e^dx. 

We recall that for a Gaussian measure with zero mean value all odd momenta 
are equal to zero. This is an important point of our considerations. This 
imply that the first nonzero contribution to the classical average is given by 
the second derivative - quadratic term. Disappearance of the third order term 
implies the asymptotics 0(a 2 ). We now estimate the remainder to obtain this 
asymptotics: 

\R(f^)\ [°° x 4 ( f\l-9fe r ^ e d9)e^dx. 

4IV27T J-oo VJ J 

Since we consider a as a small parameter, we can assume that |cr(/x)| < 1 in 
the exponential function. Thus: 



/OO 2 
x 4 e r W~^dx. 
■oo 



Since a 2 (fx) = a + 0(a 2 ), we have that R(f, /i) = 0(a 2 ), a — > 0. 

We consider the dispersion cr 2 (fi) as the intensity of fluctuations in the 
ensemble of systems. We define the relative average with respect to the inten- 
sity of fluctuations by normalizing the average by the main term - namely, 
a - in the intensity of fluctuations: 
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Of course, (f)^ is also a function of the parameter a 



</>„(„) = <!>> <»> 



a 



Corollary 4.1. Let f E V(Q) and let fj, E S%(Q). Then 

no) 



In particular, 



a->0 

Proposition 4.1. We have: 



(f)» = L -k- L + 0(a). (25) 



li"i</») = ^P-. (26) 



(f)» = —^f + 0(a). (27) 

Remark 4.1 (About 1/2) The second term in the Taylor formula gives the 
factor 1 /2 which looks rather bothering in our asymptotic formula for the classical 
average. This factor will disappear in the complex representation and the formula 
will become nicer. 

We have shown that gives the approximation of the (classical) rel- 
ative average. The precision of such an approximation is a. If the level of 
development of measurement technology is such that all contributions of the 
magnitude a are neglected in measurements, then averages can be calculated 
by using the following simple rule: 

/Aapprox [</>M 1 appr ° X / (0) / O o\ 

(/) - "b^rJ = ^r- (28) 

At the first sight such averages have nothing to do with classical averages 
given by integrals. There could be even presented an interpretation of physics 
claiming that rules of classical probability theory are violated and relating the 
exotic rule (|28)l for calculating of averages to special features of systems under 
consideration (and not to a special approximation procedure for averages). 

Finally, we remark that calculation of averages by (|2*Hj) is essentially sim- 
pler than classical probabilistic averages given by Lebesgue integrals. 
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5 Taylor approximation of classical averages: 
multidimensional case 

States are vectors q G Q = R m ; statistical states are Gaussian distributions 
with the zero mean value and the dispersion <J 2 (n) = a + 0(a 2 ). Denote this 
class of probabilities by the symbol Sq(Q). We introduce the scalar product 
and norm on Q : 

m m 
3=1 3=1 

If a Gaussian measure \i is nondegenerate (so the measure of any open set is 
positive), then 

e -W^hti dq 

^ ~ v/(2tt)™ det B' 

where B is a positive operator (we consider everywhere only Gaussian mea- 
sures with zero mean values). If ji e Sq(Q) and nondegenerate, then 

W v /(27r) m det J B > /R m " y " V ' 

In the general case the easiest way to define a Gaussian measure is to use its 
Fourier transform: 

= f e ^' q) dfi(q) = e -^ Bq ' q \ 
Jr™ 

where B = cov ji is the covariance operator: 



(56,6)= / (6,?) (Z2,q)dfi(q). 

</R m 



We remark that by definition a covariance operator is positively defined and 
symmetric. 

Lemma 5.1. Let /i be a Gaussian measure with the zero mean value and 
let A be a symmetric operator. Then 



[ (Aq,q) dfi(q) = Tr BA, (29) 

JR" 1 



where B = cov ji. 
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To prove this lemma we should just expand the quadratic form (Aq, q) 
with respect to an orthonormal basis. 

Corollary 5.1. We have 

a 2 (p) = / ||g|| 2 d/i(q) = Tr B. (30) 
Thus, for fi G S%(Q), 

Tr cov p = a + 0(a 2 ). 

We now define a class of physical variables - V(Q) : a) / G C°°(R m ); b) 
/(0) = 0;c)||/( 4 )(g)||<c / e r /NI,c / ,r / >0. 

For a function / : R m — ► R, its nth derivative is a (symmetric) n-linear 
functional, f^ n \q) '■ R m x ... x R m — > R. The norm of this functional is given 
by 

\\f (n \q)\\= sup \f^(q)(h 1: ...,h n )\. 

The norm can be estimated by partial derivatives: 

d n f(q) 



\\f {n \q)\\< max 

ai+...+a n =n 



dq" 1 ...dq% ri 



It is easy to generalize Lemma 4.1 to the multidimensional case. 

Thus we have defined the following classical statistical model: N c \ ass = 
(Sg(Q),V(Q)). 

Lemma 5.2. Let f G V(Q) and let p G S'g(Q). Then 

< / > M (a) = [ f(q) dp(q) = ^Tr pf"(0) + 0(« 2 ), (31) 
Jr" 1 * 

where p is a density operator; in fact, p = cov p/a. 

Proof. By using the scaling of the state variable (|22|) and by expanding 
f(a(p)x) on the basis of the fourth order Taylor formula with the integral 
remainder we obtain: 



<f>, (a) = ^Tr pf"(0) (32) 



^^-J (/ 1(1 " e)3 f^faWb'M'V) dd ) d ^ x ) 



15 



where /i sca i is a normalized Gaussian measure - the image of fi under the 
scaling (|22j) . We now estimate the remainder: 

w.aoi ^^rj Km ii^ii 4 ( jf c 1 - e ^ 3 eMfl)Me dd ) d ^^- 

Thus 

W,V)\ < C'a\p) [ ||x|| 4 e r W\ dp scal (x). 



R" 



We have that /i) = 0(a 2 ), a — >• 0. 

Corollary 5.2. Let f G V(Q) and let p G Sq(Q) be nondegenerate. Then 

1 f „, , icr-i-^ , a. 

A/(2vr) m det5 7 R " 

where p = B/a. 

As in the one-dimensional case, we introduce the relative average: 



/ /(?) e*< fl_1 ™> dq = ^Tr pf"(0) + 0(a 2 ). (33) 



ffl- — = TTf7?t + W. 

In the case of a nondegenerate Gaussian measure we have: 

L f(q) es( B-1 M) dg 

Corollary 5.3. Lei / G V(Q) and Zei /i G Tnen 

(/), = iTr P r(0) + O(a). (34) 

Thus if one neglects by terms of the magnitude a, it is possible to use the 
following approximative calculus of averages: 

(/)— = \^ M> (35) 

where A = /"(0) and p = cov p sca ,\- This is nothing else than the von Neu- 
mann trace formula for quantum averages, see (|13|). To proceed more for- 
mally, we consider maps: 

T : Sg(Q) - DW(m),p = T(/i) = cov /i scal ; (36) 
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T : V(Q)) - 4 r) (m)), A = T(f) = /"(O) (37) 

(we recall that Hessian is always a symmetric matrix). 

Theorem 5.1. The maps A 3 6]) . Jff7| ) project the classical statistical model 
N dass = (S%(Q),V{Q)) onto the quantum model A^Lt = (DW(m),£[ r) (m)) 
in such a way that classical and quantum averages are coupled by the asymp- 
totic equality: 

0% = \< T(f) >tm +0(a). (38) 

6 Degenerate Gaussian measures and "pure 
states" 

Consider a Gaussian measure p G Sq(Q) which is concentrated on a linear 
subspace Qq of Q. So it is nondegenerate on Qq. Here P : Q — ► Qo is the 
orthogonal projector onto Qo- Denote by So the co variance matrix of the 
restriction of p onto Qq. Then B > and i? = cov p = PB P. Let us now 
make the scaling (J22j) . Then p = cov p sca i = PpoP, where p — cov /^scailQo- 

Thus for any symmetric matrix A we have Tr pA = Tr p (PAP). Suppose 
now that p = dir ^q , where I : Qq — > Qo is the unit operator. Then 
Tr pA = dim 1 q o Tr PAP. We are especially interested in measures concentrated 
on one dimensional subspaces, Qq = Qy — {q — c^f,c G R}, where \l/ has 
the norm one. Here P = is the one dimensional projector P^ : Q ^ 
Qy, P^4> = ((p, and hence 

Tr pA = Tr P*AP* = (Atf, (39) 

Denote a probability p GG Sq(Q) which is concentrated on the one di- 
mensional subspace by the symbol p$. We obtained the following simple 
result: 

Proposition 6.1. For any p^, we have 



(f),* = l(f"m,*) + 0(a). (JO) 



Thus approximately we have: 



</u = ^*,n ^ = /"(o). (4i) 
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But the right-hand side of this equality is nothing else than the well known 
quantum formula for the average of the quantum observable A with respect 
to the pure state [3]. In our approach this quantum formula arose as the 
approximation of the classical average with respect to a Gaussian ensemble. 
The only distinguishing feature of such an ensemble is that all systems have 
states proportional to the vector \1/ (with the probability one). Of course, 
one can consider the projective space and then all those systems will have 
the same coordinate. However, real coordinates of systems are different. 

Conclusion. Quantum averages with respect so called pure states can 
be easily reproduced as approximations of ordinary ensemble averages with 
respect to one dimensional Gaussian distributions. 

7 Prequantum phase space — the two dimen- 
sional case 

In previous sections we considered the prequantum toy model in that the 
phase space structure was not taken into account. The corresponding quan- 
tum model was over the reals, see also [TT], On the other hand, physical 
reality is described by the classical phase space mechanics and the complex 
quantum mechanics. We shall see that it is possible to create a prequan- 
tum phase space model reproducing the complex quantum mechanics. The 
crucial point is that classical variables and statistical states - functions and 
measures on phase space - should be invariant with respect to a special group 
of transformations of phase space. 

This fundamental prequantum group is very simple - the special orthogonal 
group SO (2), the group of rotations of phase space. 

States of systems are now represented by points ip = (q,p) G Q = Q x P, 
where Q — P — R. Here the q is the position and the p is momentum, 
so Q denotes phase space. Statistical states are represented by Gaussian 
S'0(2)-invariant measures having zero mean value and dispersion 

a 2 (/i) = 2a + 0(a 2 ); (42) 

physical variables are by S'0(2)-invariant maps, / : f2 — > R, which satisfy 
conditions a), b), c) specifying variables in the real case. Denote these classes 
of measures and functions, respectively, Sq(Q\SO(2)) and V(Q\SO(2)). 
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The appearance of the factor 2 has the following motivation: there are two 
contributions into fluctuations - fluctuations of positions and momenta. We 
shall see that they are equally distributed. Therefore it is natural to consider 
as a small parameter of the model the dispersion of e.g. the g-fluctuations 
(which equals to the dispersion of the ^fluctuations). 

We consider the classical model iV class = (S%{Q\SO{2)),V{n\SO{2)). As 
in the real case, we can obtain the asymptotic expansion of the classical 
averages, see (|34p. However, in quantum mechanics we consider the complex 
structure. We would like to recover it in our classical model. To do this, we 
shall study in more detail properties of classical probabilities and variables. 

A measure /i is invariant if for any u G 5*0(2) : 

I f(uq)dti(q) = [ f(q)dfi(q). (43) 
Jr. 2 J-r? 

For a Gaussian measure fi with the covariance matrix B, this is equivalent 
to the condition: 

[u,B] = 0, u G SO (2). (44) 

Let / be a two times differentiate invariant map, so f{uip) = f{ip)-, for 
any u G SO (2). By representing 

/ cos 6 — sin 9 \ , A _. 

u = u = . a „ , (45) 
y sin ft cosy J 

we have that 

f(cos9q — sin9p, smOq + cosOp) = f(q,p). (46) 

This is a rather strong constraint determining a very special class of maps. In 
particular, we obtain: u*V f{uip) = V f(ip) and u* f"{uip)u = f"(ip). Hence 
ii*V/(0) = V/(0) for any rotation, and thus 

V/(0) = (47) 

and 

[f (0),H = 0, ueSO{2). (48) 

It is convenient to introduce the commutator of the set SO(2) in the 
algebra of all two by two matrices M^ r \2) : 

SO' (2) = {Ae MW(2) : [A,u] = 0,u G SO{2)} 
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We remark that a generator of 5*0(2) can be chosen as the symplectic oper- 
ator: 

' 1 
-1 



J 



Therefore the commutator of SO' (2) coincides with the commutator of J : 
{J}' = {AeM^(2):[A,J]=0}. 

Proposition 7.1. Let fi G S^(n\S0(2))and let f G V(n\S0(2)). Then 
B = cov n and A = /"(O) belong to SO' '(2). 

Lemma 7.1. A matrix A belongs to the commutator SO' (2) z/f 

s / ) • w 

If A is also symmetric, then it is diagonal: A = ^ ^ ^ ^ . In particular, 
its trace is given by 

TrA = 2R (50) 
Thus if pi G Sq(Q.\SO(2)), then its covariance matrix is diagonal B = 
q 5 ^ ' wnere 26 = a + 0(a 2 ). Fluctuations of the coordinate g and the 
momentum p are independent and equally distributed: 

dub) = irb eM ~ q -i^ }dq - 

Denote the marginal distributions of \x by the symbols /i q and /i p , respec- 
tively. Then 

1 f +co q 2 1 f +co p 2 



Hence 



(T 2 (// 3 ) = a 2 (/i p ) = ^cr 2 (fi) = a + 0(a 2 ). 



Proposition 7.2. Let f G V(Q\SO(2)). Then all its odd derivatives at 
the point qo = are equal to zero. 

Proof. For any vector G f2, we have f( 2n+1 \uip)(u<j), ...,m0) 
_ y(2n+i)^^0^ •••,0), for any rotation w. We choose u = J. Then: 

(-i) 2n+ V (2n+1) (o)(0,...,0) = ; (2 " +1) (o)(J 2 0,..., JV) 
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= /( 2 " +1 )(O)(J0, J<j>) = / (2 " +1) (O)(0, 0). 
Thus /( 2n+1 )(O)(0, 0) = for any vector (pen. 

For a function / G V(Cl\SO(2)), its Hessian has the form f"(0) = ^ ^ ^ 
where i? G R, and hence: 



2^ 



/(,,,) = *£+£> +0 ( a >). 

We remark that, in spite of the coincidence of commutators, the 50(2)- 
invariance is not equivalent to the J-invariance (the later was used as the 
basis of the theory in [T5]). 

Example 7.1. Let f(q,p) = q 3 p—qp 3 = qp(q 2 —p 2 ). Then f(Jip) — /(^)- 

But take 9 = vr/4. Here u = \ ( \ ~ l ) . Hence «(j) = ([^ ^ ) ' 

Thus f(wp) = (q — p)(q+p)(q—p — q—p)(q—p + q+p)/16 = —qp((q 2 —p 2 )/4:. 

We are now completely ready to recover the complex structure of quantum 
mechanics. By Lemma 7.1 any matrix belonging S'0 / (2) can be represented 
in the form: A = R I + S (-J). By mapping / into 1 and (—J) into % we 
obtain a map of the commutator SO'{2) onto the set of complex numbers 
C : 

j : S0'(2) ^ C,z = j(A) = R + iS. (51) 
This is the isomorphism of two fields. 

' R 



In particular, a symmetric matrix A = I ^ ^ J is represented by the 

real number j(A) = R. This is the operator of multiplication by R. The 
trace of this operator in the one dimensional complex space C (with the 
scalar product, (z,w) = zw) equals R. By (|5Uj) we obtain that 

Tr A = 2Tr j(A), (52) 

where at the left-hand side we have the real trace and at the right-hand 
side - the complex trace. Now we can write the basic asymptotic equality 
for averages in the complex form. In the funny way the Taylor factor | 
disappears through the transition form the real to complex structure, see 
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Lemma 7.2. Let f G V(Q\S0(2)) and let // G S%(Sl\S0(2)). Then 



<f >IM (a)= f(q,p)d^q,p) = aj(f"(0)) + O(a 2 ). (53) 
Jr 2 

Proof. We make the scaling of the state variable: 

* = ^ (54) 

Then the image of /i is again a Gaussian measure, say /i sca i, having the 
dispersion a 2 (/i sca i) = 2. Set £> = cov yU sca i- In the two dimensional case 

D = f q J J and Tr D = 2. We now have: 

<f>, (a) = ^Tr D f"(0) + 0(a 2 ). (55) 

Thus 

</>/.(«) = ^ J'P) J(f (0)) + 0(a 2 ). (56) 
Finally, we note that in the two dimensional case: j(D) = 1. Thus we obtain: 

</>,(«) = ^j(f(0)) + O(a 2 ), (57) 

and hence (J53j) . 

We recall that in the one dimensional quantum mechanics there is just 
one "density matrix" , namely, p = 1 G R. 

It is convenient to consider the renormalization of averages by the main 
term in the intensities of fluctuations of the coordinate and momenta: = 
< ^ >M . Then we get: 

(f)»(a)=3(f(0)) + O(a). (58) 



8 Prequantum phase space — multidimensional 

case 

States of systems are now represented by points ip = (q,p) e Q = Q x P, 
where Q = P = R m . Here the q = (qi, ...,q n ) is the position and the p = 
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(pi, ...,p n ) is momentum, so Q denotes phase space. Let us consider the 
canonical representation of the group 5*0(2) in the phase space Q = Q x P : 

( cos 9 I — sm9 I \ 
u = u e = . a Q , (59) 

y sm v 1 cos a 1 J 

where I is the unit matrix from M^(m). The corresponding group of R- 
linear operators (or 2m x 2m matrices) we denote by the symbol SO m (2). 

The classical model N dass = (S%(Q\SO m (2)),V(n\SO m (2))) in defined 
in the same way as in the two dimensional case. A Gaussian measure is 
invariant iff its covariance operator belongs to the commutator S0' m (2) = 
{A G M (r) (2m) : [A, u] = 0, u G SO m (2)}. If a smooth function / is invariant 
then all its odd derivatives equal to zero and the second derivative belong 
to the SO' m {2). A matrix A G SO' m {2) if it has the form flU, where R, S G 
M^ r \m). In contrast to the two dimensional symmetric matrix from 

SO' m {2) can be nondiagonal. It has the form (jUIJ), where R* = R and S* = 
-S. 

There is a natural map (generalizing the map j : SO' (2) — > C) of the 
commutator 5'0^ n (2) onto the set of complex matrices M^(m) : 

j : 50^(2) -> M^(m), z = = i2 + iS 1 . (60) 

This is the isomorphism of two rings. Symmetric matrices are mapped onto 
symmetric matrices. Let us denote real and complex conjugations by * and 
*, respectively. We have (R + iS)* = R* — iS* = R + iS. We also remark 

that for a symmetric complex matrix: 

Tr j(A) = Tr (R + iS) = Tr R= - Tr A. (61) 

Lemma 8.1. Let f G V(n\SO m (2)) and let p G S%(Q\SO m (2)). Then 

<f> ll (a)= a Trp j(f"(0)) + 0(a 2 ), (62) 

where p G T>( c \m). 

Proof. We make the scaling and get the /i sca i with D = cov /i sca i, 
and Tr D = 2. We set p = j(D), here Tr p = (Tr D)/2 = 1 and p G P (c) (m). 
Finally 

< / >m (a) = ^ Tr j(D) j(f"(0)) + 0(a 2 ) (63) 
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implies (JB2J) . 

We now modify the classical— > quantum projections, (J37J) . to make 
them consistent with the complex structure: 

T : S%(tt\SO m (2)) -> P (C, H, p = = j(cov /i scal ); (64) 

T : V(fl|50 m (2)) - 4 C) M, A = r(/) = j(f)(0) (65) 

Theorem 8.1. The maps < f^3)) project the classical statistical model 
iVciass = (SG( fi |S°m(2)), V(Q|SO m (2))) onto £/te quantum model = 

(V^ (m) , (m)) in such a way that classical and quantum averages are 
coupled by the asymptotic equality: 

(/)„ =< T(f) > TM +0(«). (66) 

9 Prequantum phase space 

States of systems are now represented by points ip = (q,p) e Q = Q x P, 
where Q = P = H and H is a real (separable Hilbert space) with the scalar 
product (■, ■) and the corresponding norm || • ||. Here the q G H is the position 
and the p G H is momentum, so f2 denotes phase space. The real Hilbert 
space structure on Q is given by the scalar product 

(i>i,ip2) = (91,92) + (PuP2)- (67) 

In physics H = L 2 (R 3 ) is the space of square integrable functions. Thus 
both position and momentum are functions of x G R 3 . These are simply 
classical fields. A point of such a phase space is a classical vector field 
ip(x) = (q(x),p(x)). 

Let us consider the canonical representation of the group SO (2) in the 
phase space Q = Q x P, see ()59|). where / : H — > H is the unit operator. 
The corresponding group of continuous R-linear operators we denote by the 
symbol SO H {2). 

The classical model 

iV cla ss = (S%(Q\SO H (2)), V{Q\SO H {2))) 
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in denned in the same way as in the finite-dimensional case. We just re- 
call a few basic notions from theory of differentiable functions and Gaussian 
measures on infinite-dimensional spaces. 

Let /ibea a-additive Gaussian measure on the a-field F of Borel subsets 
of fl, see [32]. This measure is determined by its covariance operator B : 
fl — > fl and mean value m = e fl. For example, B and m determines the 
Fourier transform of p : 

p( y ) = I e iM) dfji(tl>) = e ^ By > y W m > y \ yefl. 
Jn 

In what follows we restrict our considerations to Gaussian measures with 
zero mean value m — 0, where (m,y) = f Q (y,if})dp(ip) = for any y G fl. 
We recall that the covariance operator B = cov p is defined by 

(By 1 ,y 2 )= / {y 1 ,ip)(y 2 ,ip)dp(ip),y 1 ,y 2 efl, 
Jn 

and has the following properties: a). B > 0, i.e., (By,y) >0,t/Gfi; b). i? is 
a bounded self-adjoint operator, B e c). 5 is a trace-class operator 

and moreover 

Tr5= / |H| 2 d/i(# 

This is dispersion c 2 (p) of the probability p. Thus o~ 2 (fi) —TrB. 

We remark that the list of properties of the covariance operator of a 
Gaussian measure differs from the list of properties of a von Neumann density 
operator [4] only by one condition: Tr p = 1, for a density operator p. 

We can easily find the Gaussian integral of a quadratic form (by using 
expansion with respect to an orthonormal basis and using our previous results 
on the finite-dimensional Gaussian integrals): J Q (Aifj,ifj)dp(if}) = Tr BA, 
where A E C s (fl). 

The differential calculus for maps / : fl — > R does not differ so much 
from the differential calculus in the finite dimensional case, / : R" — > R. 
Instead of the norm on R n , one should use the norm on fl. We consider so 
called Frechet differentiability. Here a function / is differentiable if it can be 
represented as 

/(^o + M) = fM + f'M(^) + o(A^), 

where limn^u-vo ^Ta^ = 0- Here at each point ip the derivative f'{ip) is a 
continuous linear functional on fl; so it can be identified with the element 
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f'{ip) G fi. Then we can define the second derivative as the derivative of the 
map ip — > f'(ip) and so on. A map / is differentiable n-times iff: 

/(^o + M) = /(V>o) + /'(^o)(A^) + i/ w (Vo)(AV, AV) + ... 

+V (n) (V>o)(AVv..,AVO + 0n (AV>), 
n! 

where f^ n \ipo) is a symmetric continuous n-linear form on f2 and 

lim 1MA^)|| = 
Ha^iho ||A^|| n 

For us it is important that the second derivative f"{ipo) can be represented 
by a self-adjoint operator f"(ip )(u,v) = (f"(tp )u,v),u,v G Q. We remark 
that for tp = we have: 

/(vo = /(o) + f(o)(v) + Jr(o)(^, v) + - + -'^ + 

2 n! 

As in the finite-dimensional case the reminder can be represented in the 
integral form. 

A Gaussian measure is invariant iff its covariance operator belongs to the 
commutator 

SO' H {2) = {Ae : [A,u] = 0,m £ SO H {2)}. 

If a smooth function / : Q — > R is S'O^ (2)-invariant then all its odd deriva- 
tives equal to zero and the second derivative belong to the SO' H {2). An 
operator A G SO' H {2) if it has the form flU, where R, S E C(H). And 
A G 50^(2) n £ S (Q) if it has the form (gHJ), where = # and S* = -S. 

Let us now consider the complexification of the Hilbert space H : H c = 
H © iH. We denote the algebra of bounded C-linear operators A : H c — > i7 c 
by the symbol C(H C ). The set of self-adjoint operators (with respect to the 
complex scalar product) we denote by the symbol C S (H C ). 

There is a natural map (generalizing the map j : SO' m (2) — > C m ) of the 
commutator SO' H {2) onto the C{H C ) : 

j : SO^(2) - £(# c ), * = = R + iS. (68) 

This is the isomorphism of two rings. Self-adjoint operators (with respect to 
the real scalar product) are mapped onto self-adjoint operators (with respect 
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to the complex scalar product). We also remark that for a self-adjoint trace 
class operator A G SO' H {2) the equality (JoTj) coupling real and comlex traces 
holds. In the same way as in the finite-dimensional case we prove: 

Lemma 9.1. Let f G V(n\SO H (2)) and let /i G S%(n\SO H (2)). Then 
the asymptotic equality KUty) . where p G V(H C ), holds. 

We now consider the infinite-dimensional generalization of the classical— > 
quantum projections, (J53|) . (J55jl 

T : S£(n\SO H (2)) -> D(# c ), p = T(//) = j(cov /i scal ); (69) 

T : V(fi|,SO ff (2)) - £ s (tf c ), A = T(f) = j(/")(0) (70) 

Lemma 9.1 implies: 

Theorem 9.1. The maps FTD^ project the classical statistical model 
Ncia SS = (S%{n\SO H (2)),V{n\SO H {2))) onto the quantum model N quant (H c ) = 
(V(H C ), £ S (H C )) in such a way that classical and quantum averages are cou- 
pled by the asymptotic equality KUty) . 

We remark that the idea that the quantum averages can be coupled to 
integration with respect to the -^-function was discussed in a number of 
papers, see, e.g., [H)]~ [HI and [TH] (so called GAP-measures) as well as 
extended literature in the last paper. The main distinguishing feature of our 
approach is elaboration of technique of asymptotic expansion with respect 
to a small parameter, namely the dispersion of prequantum fluctuations. 
Comparing with ^J] we also mention that we (as well as Bach JH]- [TH] ) 
consider the linear space integration and not integration over the unit sphere. 
This is not simply a technical deviation, but it implies a totally new viewpoint 
to quantum pure states, see the next section. 

10 Gaussian measures corresponding to pure 
quantum states 

We now generalize considerations of section 6 to the complex infinite-dimensional 
case. Let \l/ be a pure quantum state: \1/ G H c , = 1. We define a Gaus- 
sian measure /i^ which is concentrated on the one dimensional complex space 
(so the real plane) 11$ = {ip G H c ■ if) = c G C} : the average of the /i$ 
is equal to zero and the complexification of its covariance operator: 
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The following facts about can be obtained through direct computations 
and Theorem 9.1: 



Proposition 10.1. For any pure quantum state we have: 

a) , G Sg(n\SO H (2)); 

b) . T(ijl 9 ) = *® 

c) . (/)„ =< j(f"{Q))% * > +0(a), f e V(Sl\SO H (2)). 



11 Hamilton-Schrodinger dynamics 

States of systems with the infinite number of degrees of freedom - classical 
fields - are represented by points ip = (q,p) G fi; evolution of a state is 
described by the Hamiltonian equations. We consider a quadratic Hamilton 
function: Tl(q,p) = |(H-0, ip), where H : Q — > Q is an arbitrary symmetric 
(bounded) operator; the Hamiltonian equations have the form: 

g = H 2 ig + H 22 p, p = -(H u q + H 12 p), 



or 



j> = (J j = mi, (7i) 

Thus quadratic Hamilton functions induce linear Hamilton equations. From 
(JZU) we get 

ijj(t) = U t ip, where U t = e jm . 
The map Utip is a linear Hamiltonian flow on the phase space Q. Let us con- 

(R T \ 
rp jl J ■ This operator 

defines the quadratic Hamilton function 

H(q,p) = ^[(Rp,p) + 2(Tp,q) + (Rq,q)], 

where the operator R is symmetric and the operator T is skew symmetric. 
Corresponding Hamiltonian equations have the form 

q = Rp-Tq, p= -{Rq + Tp). (72) 

We point out that for a 50#(2)-invariant Hamilton function, the Hamilto- 
nian flow Ut G SO' H (2). By considering the complex structure on the infinite- 
dimensional phase space Q we write the Hamiltonian equations (|71j) in the 
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form of the Schodinger equation on H c : 

(73) 

Its solution has the following complex representation: tp(t) = U t ip, U t = 
e _lHt . We consider the Planck system of units in that h — 1. This is the 
complex representation of flows corresponding to quadratic SO h {2) -invariant 
Hamilton functions. 

By choosing H = L2(R n ) we see that the interpretation of the solution of 
this equation coincides with the original interpretation of Schrodinger - this 
is a classical field 

4>(t,x) = (q(t,x),p(t,x). 
Example 11.1. Let us consider an important class of Hamilton functions 

n(q,p) = ^[(Rp,p) + (Rq,q)], (74) 

where R is a symmetric operator. The corresponding Hamiltonian equations 
have the form: 

q = Rp, p = —Rq. (75) 

We now choose H = L 2 (R 3 ), so q(x) and p(x) are components of the vector- 
field ip(x) = (q(x) , p(x)) . We can call fields q(x) and p(x) mutually inducing. 
The field p(x) induces dynamics of the field q(x) and vice versa, cf. with 
electric and magnetic components, q(x) = E(x) and p(x) = B(x), of the 
electromagnetic field, cf. Einstein and Infeld [20]; P- 148: "Every change 
of an electric field produces a magnetic field; every change of this magnetic field 
produces an electric field; every change of and so on." We can write the form 
(JZH) as 

ft(?,p) = o/ R(x,y)[q(x)q(y) + p(x)p(y)}dxdy (76) 

1 </R6 

or 

^C0) = o / R{x,y)^{x)'4>{y)dxdy, (77) 

where R(x, y) = R(y, x) is in general a distribution on R 6 . We call such a 
kernel R(x,y) a self-interaction potential for the background field ip(x) = 
(q(x),p(x)) We point out that R(x, y) induces a self-interaction of each com- 
ponent of the i[}(x), but there is no cross-interaction between components 
q(x) and p(x) of the vector-field ip{x). 
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12 Invariant Gaussian measures of the Hamilton- 



Schrodinger dynamics and stationary pure 
states 



All Gaussian measures considered in this section are supposed to be SOh(2)- 
invariant. As we have seen, In our approach so called pure states \&, 1 1^| | = 1, 
are labels for Gaussian measures concentrated on one dimensional (complex) 
subspaces of the infinite-dimensional phase-space Q. In this section we 
study the case of so called stationary (pure) states in more detail. The en- 
sealing does not play any role in present considerations. Therefore we shall 
not take it into account. We consider a pure state | | = 1, as the label for 
the Gaussian measure having the zero mean value and the complexification 
of the covariance operator 



Theorem 12.1. Let v be a Gaussian measure (with zero mean value) 
concentrated on the one- dimensional (complex) subspace corresponding to a 
normalized vector ^ . Then v is invariant with respect to the unitary dynamics 
U t = e~ ttH , where H : Q — > Q is a bounded self-adjoint operator, iff ^ is an 
eigenvector of H. 

Proof. A). Let H\I> = A^. The Gaussian measure U*u has the complex- 
ification of the covariance operator 

j(p t ) = U t (^ <g> #)[/ t * = U t V <g> U t V = e~ itx ^ <g> e~ itx ^ = <g> 

Since all measures under consideration are Gaussian, this implies that U^u = 
v. Thus v is an invariant measure. 

B). Let Ufv = v and v = for some ||^|| = 1. We have that U t ^ <S> 
Ut^ — \& <S> ^. Thus, for any ipi, tp 2 € we have 



Let us set ^2 = ^ - We obtain: c(t)U t ^) = (V'i,*), where c(t) = 
(U t V,V). Thus c(t)Ut$ = fr.We point out that c(0) = ||*|| 2 = 1. Thus 
d(0)^-iH^ = 0, or = -ic'(O)*. Thus * is an eigenve ctor o f H with the 
eigenvalue -i c'(O). We remark that c'(0) = -i(Hty, so c'(0) = i(H\P, 
Hence, A = -ic'(O) = (H, , 




(78) 
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Conclusion. In PCSFT stationary states of the quantum Hamiltonian 
(represented by a bounded self-adjoint operator H) are labels for Gaussian 
one- dimensional measures (with the zero mean value) that are invariant with 
respect to the Schrddinger dynamics Ut = e~ tm . 

We now describe all possible Gaussian measures which are Ut- invariant. 

Theorem 12.2. Let H be a bounded self-adjoint operator with purely 
discrete nondegenerate spectrum: H^ k = X k ^ k , so {^ k } is an orthonormal 
basis consisting of eigenvectors o/H. Then any U t -invariant Gaussian mea- 
sure v ( with the zero mean value ) has the complexification of the covariance 
operator: 

oo 

Hp) = ^c fc ^ fc ® ^ k ,c k > 0, (79) 

k=l 

and vice versa. 

Proof. A). Let j(p) has the form (J7HD. Then 

oo 

j(pu?») = U t j(p)U; = Y,Cke- lX ^ k ®e- iX ^ k =j(p). (80) 

k=l 

Since measures are Gaussian, this implies that U*v — v for any t. 

B). Let U^v — v for any t. We remark that the complexification of any 
covariance operator p can be represented in the form: 

oo 

j(p) = mo'(p)*fe» ® + 2o'(p)**> ® *«■ ( 8i ) 

k=l k^l 

We shall show that (j(p)^> k ,^/j) = for k ^ j. Denote the operator corre- 
sponding to Ylikj^j by Z. We have 

{UtZUt^M = Y,^P)^k^ m )e lt ^ Xk \^ k ^2)^i^ m ) = {ZMi). 

k^m 

(82) 

Set = ^j,ip2 = ^k- Then 

(U t ZU^ m , * fc ) = {j{p)* k , V m )e lt ^-^ = (j(p)* k , * m >. (83) 
Thus 0'(p)*fc,*m> = 0,fc^m. 
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13 Stability of hydrogen atom 



As we have seen, in PCSFT so called stationary (pure) states of quantum me- 
chanics can play the role of labels for Gaussian measures (which are SOh(2)- 
invariant and have zero mean value) that are U t - invariant. We now apply 
our standard a-scaling argument and we see that a stationary state \1/ is 
a label for the Gaussian measure ^ with j{p^) = ct^ <8> VP. This mea- 
sure is concentrated on one-dimensional (complex) subspace LT^ of phase 
space VL. Therefore each realization of an element of the Gaussian ensemble 
of classical fields corresponding to the statistical state gives us the field 
of the shape ^(x), but magnitudes of these fields vary from one realization 
to another. But by the well known Chebyshov inequality probability that 
£(^) = J R 3 |\l/(a;)| 2 efo; is large is negligibly small. 

Thus in the stationary state we have Gaussian fluctuations of very small 
magnitudes of the same shape ty(x). In PCSFT a stationary quantum state 
can not be identified with a stationary classical field, but only with an ensem- 
ble of fields having the same shape ty(x). Let us now compare descriptions 
of dynamics of electron in hydrogen atom given by quantum mechanics and 
our prequantum field theory. 

In quantum mechanics stationary bound states of hydrogen atom are of 
the form: 

V nlm (r, 9, 0) = c n , l R l L 2 ^(R)e' R / 2 Y l m (9, 0), (84) 

where R — and a = is a characteristic length for the atom (Bohr 

radius). We are mainly interested in the presence of the component e~ R / 2 . 

In PCSFT this stationary bound state is nothing else, but the label for 
the Gaussian measure , which is concentrated on the subspace Q*, , . 
Thus PCSFT says that "electron in atom" is nothing else than Gaussian 
fluctuations of a certain classical field, namely the field ^f n i m (r, 9, 0) : 

Vw(r, 9, 0; V) = lW^nim(r, 9, 0), (85) 

where ^(ip) is the C- valued Gaussian random variable: E^y = 0, E\^\ 2 = a. 

The intensity of the field ^ n i m (r, 6, 0, ip) varies, but the shape is the same. 
Therefore this random field does not produce any significant effect for large 
R (since e~ R / 2 eliminates such effects). 

Thus in PCSFT the hydrogen atom stable, since the prequantum random 
fields ipnimij'i 0, 0; -0) have a special shape (decreasing exponentially R — > oo). 
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14 Appendixes 



14.1 Classical representation for spin operators 



The Pauli matrices are a set of 2 x 2 complex Hermitian and unitary matrices. 



They are: o\ = 



1 

1 



, <r 2 = 



. Let H c = C 2 with 



—i\ 1 

t o ) ' as ~ [o -1 
the complex coordinates z — (z±, z 2 ), Zj = qj + ipj,j — 1, 2, and = R 2 xR ; 
with the real coordinates uj = (qi,q2,Pi,P2)- We consider spin operators 

= E? a j°j : c 2 

^ r) = r>i) = 



C 2 , a = (01,02,03). Let us consider real matrices 



a 



(r) 



0-1 
01 



/o 


1 





o\ 


1 




















1 


\o 





1 


0/ 


(r) 




'0-3 





3 


■( 





0-3 



(r) 
ff 2 



—%a 2 
%o 2 



/o 
















1 








1 








V- 1 








0/ 



/1 








o\ 





-1 














1 













-v 



We remark that these are not Dirac matrices. We set 

3 



(7W(a)=rXa))=E°^. 



(r) 



R 



R 4 , a = (ai, a 2 , a 3 ). 



±(aW(a)u,u). Then 
R 2 x 



and consider classical random spin-variables: f a (oo) 
T(fa) = cr(a) and for any S^^-invariant Gaussian measure /ionfi 
R 2 with dispersion a+0(a) we have: ^ J R4 f a (uj)d^(uo) = Tr j(p)(j(a) + 0(a), 
where p is the covariance operator of ^a-scaling of the Gaussian measure p. 
For example, 



1 



(qiq2+pm)dli(qi,q 2 ,Pi,P2) = Tr j(p)<7i + 0(a), 



- / (Pi52 -P2qi)dn(q!,q 2 ,pi,p 2 ) = Tr j(p)tr 2 + 0(a), 
/ (?i - 92 + Pi -p1)<M9i»92,Pi,P2) = Tr j(p)cr 3 + 0(a). 

JR 4 
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We also have: ± J R4 (qf + q\ + pj + p 2 2 )dfx(q u qi,Pi,Pi) = Tr j(p)I + 0(a) = 
T" j(p) + 0(a). By introducing vectors u>i = (qi,Pi) and uj 2 = (<Z 2 ,p 2 ) we 
rewrite these asymptotic equalities in shorter way: 

- / {uji,uj 2 )d^(uji,uj 2 ) = Tr j(p)(7i + 0(a), 

- / (Jui,u 2 )dfj,(ui, lo 2 ) = Tr j(p)<7 2 + 0(a), 
where J is the symplectic operator and, finally, 

- / (|M| 2 - ||w 2 || 2 )dAi(wi,w 2 ) = Tr j(p)cr 3 + 0(a). 
a Jr4 

Let us now consider Gaussian measures on Q = R 4 corresponding to 
pure quantum states. These are singular Gaussian measures which are con- 
centrated on 5'0 2 (2)-invariant planes in R 4 . To determine such a measure, 
we should find its covariation operator. 

Proposition 14.1. Let ^ = u + iv,u = (u±, u 2 ) e R 2 , v = (v 1, v 2 ) G R 2 

be a pure quantum state and let = \& x ^. Then T _1 (p^,) = /i^, where the 
Gaussian measure ji^, has the covariation operator By = aD^ for 

12 



( IN 
(91,92) 


\(Jgi,92) 



(91,92) 
W92W 2 
(gi,Jg2) 






(gi, Jg2) 
hif 
(91,92) 



(Jgi,g2)\ 


(01,02) 
II02II 2 ) 



Here g\ = (u±,Vi) and g 2 = (u 2 ,v 2 ) are variables which are conjugate to 
vi = (qi,Pi) andu 2 = (? 2 ,p 2 ). 

Proof. The real space realization of is given by the operator: 



u ® u + v <g) v v0u — u®v 
u®v—v®u u®u+v®v 



r\ P *) = 

We have in the chosen system of coordinates on the phase space: 

UiU 2 \ 



u x u 2 



Un 



J 



V (X) V = 



V\V 2 



y VlV 2 
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v <S> u 



UiVi U 2 V\ 

u 1 v 2 u 2 v 2 



V\Ui V 2 U\ 
v x u 2 V 2 U 2 



Hence: 



U (g) U + V (g) V 



V ® u — U ® V 



U\U 2 + v x v 2 

UiU 2 + V\V 2 «2 + V 2 

^2^1 - V 2 Ui 

u\v 2 - v\u 2 



To illustrate better correspondence between real and complex state spaces, 
we now show directly that j(Dq,) = \I/ ® \l/ for given by this Proposition. 
We have 



j(D 9 ) = 



\9i\ 



(91,92) 



(J 91,92) 

(gi,Jg2) 



^(91,92) \\92W 2 

This operator acts to a complex vector z — (zi,z 2 ) in the following way: 
A = U( D *) z )i = \\9i\\ 2 zi + [(91,92) +K J 9i,92)]z2, 

4 = CK- *)^ = [(fi-i,^) + 1(^1,^2)]^ + ll^ll 2 ^- 

On the other hand, * ® #(V) =< z, * > * = (zi^i + 22*2)*- Here 
4 = (mi - Wi)(mi + iv 1 )z 1 + (u 2 - iv 2 ){u 1 + Wi)z 2 , 
4 — ( u i ~ ivi)(u 2 + ^2)^1 + (u 2 - iv 2 )(u 2 + iv 2 )z 2 . 



Thus 



4 = ( M l + v l) z l + [(^1^2 + ViV 2 ) + i(ti2^1 - Ml^)]^, 
4 = ( M 2 + ^2)^2 + [iUiU 2 + + i(UiV 2 - W 2 ^l)]^2- 

Let us consider Gaussian measures corresponding to pure states for spin 



up and spin down, |1 >= 



and 10 >= 



Ui 



1,^2 



Vl 



v 2 



0. Thus £)<!> 



transform of the measure /i|i> is given by: 

£|i>(fi, Z2,r)i,m) = e'^l+vl) . 



. For vector |1 >, we have: 

/l 0\ 

10 

\o 0/ 



and the Fourier 



35 



This is the standard Gaussian measure on the plane q 2 = 0,p 2 = hav- 
ing the density: rf/i|i>(?i,Pi) = ~2swi+ p i). j n the same wa y ^| 0> j s the 
standard Gaussian measure on the plane qi — 0,pi — having the density: 



2-7TQ 



e 2 Q ('?2+P2). Let us now consider the Gaussian measure 



corresponding to superposition of spin up and spin down states: 



-^(|0>+e* e |l>). 



Here u\ = cos^fi = sin 9. Hence 



( cos 2 9 cos 9 sin 9 

cos # sin 9 sin 2 # 



\ 






cos 9 sin 6* 
cos 9 sin # sin 2 J 






cos 2 6 1 



and the Fourier transform of /i^ fl is given by 



(£1,6, ^1,^2) 



e 2 



[(cos 6»£i +sin 6»£ 2 ) 2 + (cos 771 +sin 6772 ) 2 ] 



Thus pure states correspond to the standard Gaussian measures concen- 
trated on planes obtained by rotations. 



14.2 Comparing with no-go theorems of von Neumann, 
Cohen-Specker and Bell 

There are no-go theorems for mathematical attempts to have a map from 
classical variables to quantum operators which preserves statistics, e.g., the- 
orems of von Neumann, Cohen-Specker and Bell, see [I], — - The no-go 
theorems say: No such map exists. In this paper we constructed such a map. 
What goes? 

Our construction does not contradict to known no-go theorems, since our 
map T does not satisfy some conditions of those theorems. An important 
condition in all such theorems is that the range of values of a classical variable 
/ should coincide with the spectrum of the corresponding quantum operator 
T(f) - "the range of values postulate." This postulate is violated in our 
framework. As we have seen, the classical spin variables are continuous and 
the quantum spin operators have discrete spectrum. Nevertheless, classical 
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averages can be approximated by quantum. Our prequantum classical sta- 
tistical model is not about observations, but about ontic reality (reality as it 
is when nobody looks at it). 

Henry Stapp pointed out [21]: "The problem, basically, is that to apply 
quantum theory, one must divide the fundamentally undefined physical world 
into two idealized parts, the observed and observing system, but the theory 
gives no adequate description of connection between these two parts. The 
probability function is a function of degrees of freedom of the microscopic 
observed system, whereas the probabilities it defines are probabilities of re- 
sponses of macroscopic measuring devices, and these responses are described 
in terms of quite different degrees of freedom." Since we do know yet from 
physics so much about features of classical — > quantum correspondence map 
T, we have the freedom to change some conditions which were postulated in 
the known no-go theorems - for example, the range of values condition. Re- 
jection of this assumption is quite natural, since, as was pointed by Stapp, a 
classical variable / and its quantum counterpart T(f) depend on completely 
different degrees of freedom. 

14.3 Is prequantum classical statistical field theory non- 
local? 

As we have seen, PCSFT does not contradict to the known no-go theorems, 
in particular, to Bell's theorem. Therefore this theory might be local. How- 
ever, it is not easy to formulate the problem of locality/no nlocality in the 
PCSTF-framework. It is not about observations. Thus we could not apply 
Bell's approach [21] to locality as locality of observations. On the other 
hand, on the ontic level PCSTF operates not with particles, but with fields. 
At the first sight, such a theory is nonlocal by its definition, since fields are 
not localized. But in field theory there was established a different viewpoint 
to locality and we know that both classical and quantum field theories are 
local. To formulate the problem of locality for PCSTF, we should proceed 
in the same way. Therefore we should develop a relativistic version of PC- 
STF. There are some technical and even ideological problems. As we know, 
relativistic quantum mechanics is not a well established theory (at least this 
is a rather common opinion). Thus it is meaningless to develop a relativistic 
variant of PCSTF which would reproduce relativistic quantum mechanics. 
The most natural way of development is to construct a kind of PCSTF not 
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for quantum mechanics, but for quantum field theory and study the problem 
of locality in such a framework. It is an interesting and complicated problem 
which will be studied in coming papers of the author. 
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